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There has been a significant interest on plasmonics in a metallic structure with very narrow gaps for studies 
of nanophotonics. However, little attention has been paid to the behavior of surface plasmons (SPs) in 
quasi-continuous metallic structures. This study observes and analyzes intermediate characteristics between 
propagating SPs (PSPs) and localized SPs (LSPs) in a quasi-continuous metallic monolayer of core-shell 
nanocubes. We reveal that, in a very narrow region of few-nanometer gaps among the nanocubes, the 
intrinsic energy bands of PSPs and LSPs intersect each other to generate two hybrid bands and an 
anti-crossing. Using a self-assembly method instead of the lithographic techniques which have several 
limitations as of now, we materialize the quasi-continuous metallic layer with plenty of nano-gaps that 
exhibit intermediate plasmonic characteristics. The intermediate plasmonic characteristics observed in this 
study will lead to interesting subjects, such as band engineering and slow SPs, in nanophotonics. 

I n macroscopic electromagnetism, a material with negligible discontinuities can be considered as a continuous 
I medium if the length scales of the discontinuities are far below the wavelength. In the surface plasmons (SPs), 
I defined as the collective oscillations of free electrons stimulated by light at an interface between metal and 
dielectric, however, the field can be significantly influenced by very narrow discontinuities which are as small as 
one ten-thousandth of the wavelength 13 . Therefore, in plasmonics, more stringent criteria are required for a 
continuous medium. To be more specific, it is expected that the continuity of a metallic structure has a relation 
with the dominant state of SPs in that structure. As is well-known, there are two independent states in surface 
plasmons (SPs): propagating SPs (PSPs) and localized SPs (LSPs). PSPs propagate along the smooth and con- 
tinuous surface of a metal layer 4 6 , while LSPs reside around a small or distinct metallic nanostructure 7-10 . 

Based on this perspective, our attention is naturally focused on the discontinuity- related intermediate char- 
acteristics of SPs. By intuition, it is expected that a transition of SPs from LSPs to PSPs will be observed if the sizes 
of the discontinuities in a metallic layer are reduced from significant to negligible lengths and eventually to zero, 
and vice versa. Since the light fields associated with SPs can be confined in a small inclusion of discontinuity of a 
metallic structure, such transition between the PSPs and LSPs may occur for very narrow discontinuities, of which 
the size is far below the diffraction limit. Although there have been several studies about the coupling between 
LSPs and PSPs, including the periodic metal particles coupled to a separated metal film 1113 or metallic holes 14-16 , 
they did not consider the transition behavior of plasmonic states in an almost-continuous metallic layer. In 
respect of the almost-continuous structure, some research groups have introduced a metal-insulator transition 17 
and tunable plasmonic property 18 at close-packed metal nanoparticles. While the previous researches broadly 
suggest the potential of a bottom-up strategy and their properties, this study can provide a detailed map of unique 
dispersion bands of SPs in one structure and concretely demonstrate their intermediate behaviors depending on 
structural parameters, as well as an in-depth understanding of the boundaries between the states of the SPs. 

This study intends to investigate the states of SPs in a metallic layer with very narrow discontinuities. The sizes 
of the discontinuities considered in this study are ranged from zero to a few nanometers. As one of the feasible 
implementations for the said structure, we propose a monolayer of metal nanocubes with thin dielectric shells. 
Also, a bottom-up approach is used, instead of the top-down process, such as lithography which is extremely 
difficult to apply, so as to build a self-assembled metallic monolayer with nano-sized discontinuities. Nano- 
discontinuities can be materialized with few-nanometer-thin dielectric shells among the neatly spaced metal 
nanocubes. Hereafter, we refer to the monolayer of nanocubes in this study as "a quasi-continuous metallic layer 
(QCML)". From the simulations and experiments of reflectometry, it is found that two distinct energy bands 
originating from the PSPs and LSPs can be observed with a phenomenon of anti-crossing in QCML. Each of the 
two bands exhibits intermediate natures of PSPs and LSPs due to their strong coupling. The dispersive nature of 
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QCML reveals the evolution or transition mechanism of the LSPs in 
the isolated metallic structure into PSPs in the continuous structure. 
In addition, the QCML provides a novel way of tailoring the disper- 
sion characteristics of the SPs or generating slow SPs by adjusting the 
size of discontinuities. 

Results 

A concept of a monolayer of core-shell nanocubes for a quasi- 
continuous metallic layer. Let us consider a monolayer of tightly- 
spaced metal nanocubes as shown in Fig. 1(a). Throughout this 
study, we use the silver nanocubes (AgNC) with silica (Si0 2 ) shells 
(AgNC@Si0 2 s) as the elements of QCML. If the shell-thickness of 
cubes is zero, the structure can be considered as a continuous metal 
layer so that PSPs can be excited via the Kretschmann configura- 
tion 46 . If the cubes are largely separated with thick shells, the layer 
can be assumed as a monolayer of particle array so that only the 
isolated LSPs can be excited 710 . Meanwhile, when the distances 
among the metal cubes are in a range from zero to the distances 
enough to isolate the cubes, intermediate plasmonic characteristics 
can be found between PSPs and LSPs due to their strong coupling in 
the QCML. For the structures with different inter-particle distances, 
we examined the dispersive characteristics of SPs via angle-resolved 
spectral reflectometry, as shown in Fig. 1 (b). Since the incident angles 
under consideration are in the total internal reflection (TIR) regime, 
the minima in the reflection curve are corresponding to the 
excitation of SP modes 4-6 . 

Intermediate plasmonic characteristics and anti-crossing. In this 
study, the results were produced using a Fourier-modal analysis 
based numerical method 19 22 . In order to make the process of 
calculation simple and efficient, we performed two-dimensional 
(2D) simulations instead of three-dimensional (3D). For the 
AgNC@Si0 2 monolayer consisting of 35 nm AgNC core and 2 nm 
Si0 2 shell, by varying the incident wavelengths from 400 nm to 
800 nm under the incident angles from 40° to 80°, the dispersion 
characteristic (0-X R map) was obtained as shown in Fig. 1(c). 

Interestingly, while only one intrinsic band for PSPs or LSPs was 
observed in a continuous case and an isolated case (insets of 
Fig. 1(a)), the two bands of SPs appeared when the shell thickness 
was in the intermediate range between the two extreme cases. In the 
conventional energy- momentum (E-k) diagram in Fig. 1(d) that is 
clockwisely rotated by 90° from the 0-A R map in Fig. 1(c), the extra- 
polated bands (the red and blue dotted lines) of PSPs and LSPs 
accord with the traditional P- and L-bands in Figs. 1(e) and (f). 
Since PSPs can be excited via the Kretschmann configuration using 
the prism, the L-band can be particularly separated from the two 
bands in Fig. 1(c) by the removal of the prism under the QCML, as 
shown in Fig. 1 (f). Further details about P- and L-bands are explained 
in Fig. SI of the Supplementary Information. In this QCML, there- 
fore, the two distinct plasmonic modes of PSPs and LSPs coexist in 
the same structure, and P- and L-bands cross each other. 
Interestingly enough, at the supposed crossing point of these two 
bands, an anti-crossing phenomenon, which is a manifestation of 
strong coupling of two modes, is observed. 

To obtain further physical insights for the intermediate character- 
istics and the anti-crossing, we examined R^-field profiles at six 
distinguishing points as shown in Fig. 1(c). In Figs. 1(g)— (I), the 
TM-polarized planewave of light is obliquely incident onto the 
metallic layer near z = 0 from the downside. In this figure, 
the supposed crossing point of two intrinsic bands is depicted as 
'O' in Figs. 1(c) and (d). When the excitation conditions are far from 
the 'O' point, the field profiles of the SPs in the intermediate structure 
resemble those of the intrinsic PSP (Figs. 3(g) and (1)) and LSP 
(Figs. 3(i) and (j)) modes. However, near the 'O' point (Figs. 3(h) 
and (k)), the field profiles for the upper (left in Fig. 1(c)) and the 
lower (right in Fig. 1(c)) energy bands in Fig. 1(d) exhibit hybridized 



modes of PSPs and LSPs. In Fig. 1(h), the fields for the LSPs in the 
metallic layer are in anti-phase with those for the PSPs on the surface. 
On the contrary, PSPs and LSPs are in-phase in Fig. 3(k). As a 
consequence, the anti-phase hybrid modes of PSP and LSP are in 
the higher (left in Fig. 1 (c)) energy band while the in-phase modes are 
in the lower (right in Fig. 1(c)) energy band. Such hybridizations of 
modes and corresponding energy relations are similar to those that 
can be found in the hybrid modes of coupled waveguides 23-26 . 

In brief, PSPs and LSPs coexisted in our QCML and anti-crossing 
was observed as an evidence of the strong coupling of these two types 
of SPs, which indicates that the PSPs and LSPs consist of the eigen- 
states in our plasmonic structure. In addition, in this study, we chose 
water (refractive index 1.33) as the surrounding material because the 
behavior of the bands is well observed under that condition. The 0-X R 
maps of the structure with air as the surrounding material are shown 
in Fig. S2 of the Supplementary Information. 

Effects of core size and shell thickness. It is important to investigate 
the intermediate range and the behavior of the two hybrid bands 
according to the structural parameters including Si0 2 shell 
thickness and AgNC core size. First, in Fig. 2(a), the increment of 
shell thickness induces a blue-shift in the location (energy level) of 
the L-band. It is because the inter-particle distance which is 
determined by shell thickness decides the resonance wavelength of 
the LSPs 27-29 . In this case (35 nm of AgNC core), the intermediate 
characteristics by the strong coupling between PSPs and LSPs can be 
observed in a very narrow region where the shell thickness is below 
8 nm, because P- and L-bands start to overlap at that thickness. 

The behavior of the bands is also affected by the size of the metal 
core as in Fig. 2(b): the core size determines mainly the relative 
strength of the PSP and LSP modes in their coupling. As the core 
size gets larger, the excitation of the PSPs gets weaker since the 
launching of PSPs with the Kretschmann configuration through a 
thick metal becomes harder. Therefore, the depth of the reflection 
dip for the PSP mode (P-band) in Fig. 2(b) gets weaker with the 
increment of the core size. On the contrary, the excitation of the 
LSPs gets weaker with the decrement of the core size. In Fig. 2(b), 
the L-band abates rapidly as the core size gets smaller because the 
scattering and absorption intensities of particles generally are pro- 
portional to the six and third power of a particle size, respectively 30,31 . 

In brief, as the sizes of dielectric shell and metal core increase, LSP 
mode becomes more dominant than PSP mode in the intermediate 
state, and vice versa. The behavior of each constituent part of two 
hybrid bands properly follows each fundamental principle of PSP 
and LSP, and that is an obvious evidence for the strong coupling 
between PSPs and LSPs. In addition, it is noteworthy that the two 
hybrid bands can be readily designed by engineering the geometrical 
parameters. 

Experimental 0-l R map of a self-assembled monolayer of AgNC@ 
Si0 2 s. According to the simulation results, in order to ascertain 
experimentally the existence of the two bands in a 0-). R map, metal 
nanocubes should have a few nanometers of narrow gaps among 
them. However, carving straight lines with a few nanometers in 
width on the metal layer is almost infeasible because the minimum 
line width of current lithography techniques is approximately 
20 nm 32,33 . Alternatively, we constructed the proposed structure 
with metallic nanocubes because the bottom-up method can 
manufacture materials in a few nanometer scales 34-37 . Cubic 
structures have a geometric advantage in forming a monolayer: the 
cubes can contact with each other or with the substrate face-to-face 
so that they can be assembled denser than other shaped structures 
even if they are randomly scattered on the substrate. They were 
synthesized with the polyol method 38,39 and enclosed with Si0 2 
shell in the Stober process 40,41 . A monolayer of AgNC@Si0 2 s was 
prepared by the self-assembled monolayer (SAM) method, which 
uses differences between solubility and surface tension of water 
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Figure 1 | Intermediate plasmonic state in QCML. ( a) Concept structures for the strong coupling between PSP and LSP in a monolayer of metallic cubes. 
Insets show the traditional 6-1r maps of PSP and LSP, respectively, (b) Simulation condition. One period includes a core-shell cube. P, d, f, 8, and R are 
period, core diameter, shell thickness, incident angle, and reflectance, respectively, (c) 6-Xr map of AgNC (35 nm) with shell (2 nm). (d) E-k diagram 
expressed by 90° rotation of the 8-a r map. (e) P-band. (f) L-band. H y -field profiles at 6 and X R of (g) 75° and 482 nm, (h) 62° and 615 nm, (i) 50° and 
637 nm, (j) 75° and 645 nm, (k) 66.5° and 660 nm, and (1) 62.3° and 750 nm. They show (g)-(i) transition from PSP to LSP and (j)-(l) transition 
from LSP to PSP. 



and organic solvent where AgNC@Si0 2 s are dispersed 42 44 . As shown 
in Figs. 3(a) and (b), AgNC@Si0 2 has approximately 43 nm of core 
size and 3.5 nm of shell thickness in the high resolution transmission 
electron microscope (HR-TEM) image and they are settled into a 
monolayer with high population density in the field emission scan- 
ning electron microscope (FE-SEM) image. AgNC orientations were 
measured from the FE-SEM image because nanocube orientations, 



including face-to-face, face-to-edge, and edge-to-edge orientations, 
have an effect on the plasmonic property 45 . As a result, we can 
confirm that more than 75% of the nanocubes contact face-to-face 
with each other, as shown in Fig. S3 of the Supplementary Infor- 
mation, and that the AgNC orientations of the experimental samples 
may be significantly in accord with the face-to-face orientation of the 
simulation. Figure 3(c) shows an illustration of the microscope setup 
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Figure 2 | Discontinuity-dependent behavior of two hybrid bands. Effect of (a) Si0 2 shell thickness and (b) AgNC core size on two hybrid bands of 9-X R 
map. 



in order to measure an incident angle-dependent reflectance spec- 
trum. The spectral reflectance of the sample was obtained with a 
spectrometer at several incident angles using a white light source 
in the Kretschmann configuration. 

Three representative cases of experimental 0-a r maps are shown 
in Fig. 4 with the corresponding simulation 0-X R maps. The three 



samples are a continuous Ag layer (thickness 43 nm) evaporated on 
glass in Figs. 4(a) and (d), AgNC@Si0 2 SAM (core size 43 nm, shell 
thickness 3.5 nm) on glass in Figs. 4(b) and (e), and AgNC@Si0 2 
randomly scattered on glass with various broad gaps in Figs. 4(c) and 
(f), respectively. They were measured for water (refractive index 
1.33) as surrounding medium and the corresponding results for air 
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Figure 3 | Experimental setup for SP excitation in AgNC@Si0 2 SAM. (a) TEM image of AgNC@Si0 2 . (b) SEM image of AgNC@Si0 2 SAM. 

(c) Schematic diagram of optical setup for measurement of incident angle-dependent reflectance spectra. Scale bars are (a) 20 nm and (b) 200 nm. 
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(refractive index 1.00) are shown in Fig. S4 of the Supplementary 
Information. The spectral reflectance data are accumulated into an 
experimental 0-a r map by scanning the incident angles between 60° 
and 75° at an interval of 1°. A series of the minimum reflectance dips 
in the experimental graphs can be regarded as a quasi-band since 
they shift with a specific tendency in wavelength. They are marked 
with round black dots on the corresponding positions in the simu- 
lation 0-A R maps in Figs. 4(d)-(f) to compare with each other. The 
single quasi-bands almost accord with P-band due to PSPs, and with 
L-band due to LSPs in Figs. 4(d) and (f), respectively. On the con- 
trary, two dips appear in one spectrum and two quasi-bands emerge 
on the experimental 0-l R map in Fig. 4(b). As the incident angle 
increases, one quasi-band near 460 nm moves toward short wave- 
length while the other quasi-band near 600 nm keeps its position. 
The experimental 0-?. R map in Fig. 4(b) can be the evidence for the 
existence of the strong coupling between PSPs and LSPs because the 
number and behavior of the quasi-bands are similar to the calculated 
6-l R map in Fig. 4(e). 

Statistical feedback of spatial distribution of AgNC@Si0 2 s to 
simulation. We can reduce the differences between the experimental 



and calculated 0-k R maps by feeding back statistical data for a gap 
distribution of AgNC@Si0 2 SAM to the prior simulation. It cannot 
be helped that a SAM of nanocubes has a certain distribution of gaps 
between them in a real sample while the nanocubes can be set at 
exact lattice points with the same distance of gap in a simulation. The 
limitation of the self-assembly method brings about differences in 
distance among the cubes: for example, two bands of Fig. 4(b) are 
placed at a shorter wavelength than those of Fig. 4(e) and the slope of 
the left band is different between Figs. 4(b) and (e). In order to solve 
the problem, we measured the inter-particle distances from the FE- 
SEM image of Fig. 3(b), as shown in Fig. 5(a). The average inter- 
particle distance was 49.3 nm and the standard deviation was 7.47 
nm. One hundred random numbers were generated in a normal 
distribution with the average and standard deviation as a distri- 
bution chart shown in Fig. 5(b), and then inputted in one period 
as the schematic diagram in Fig. 5(c). Namely, the calculation was 
done for one hundred nanocubes randomly and normally distri- 
buted in one period. 

The recalculated 0-). R map is demonstrated in Fig. 5(f) along with 
the initially calculated map in Fig. 5(d) and the experimental map in 
Fig. 5(e) for comparison. The corresponding results for air (refractive 
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Figure 4 | Experimental and simulation 0-k R maps, (a)-(c) experimental and (d)-(f) calculated 6-k R maps for AgNC (43 nm) with shell thickness 
of (a) and (d) 0 nm, (b) and (e) 3.5 nm, and (c) and (f) 21.5 nm. The positions of the dips on the experimental maps of (a)-(c) are marked with round 
black dots at the same positions on the simulation maps of (d)-(f). The whole regions in the experimental maps are denoted by dotted boxes on the 
simulation maps. 
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index 1.00) as surrounding medium are shown in Fig. S5 of the 
Supplementary Information. When the spatial distribution of the 
particles is applied to the simulation condition, two bands become 
broad and the reflectance increases slightly while the overall tend- 
ency of 0-l R map remains the same. In particular, the blue-shift of the 
left band is bigger than that of the right band because the short 
wavelength of light can identify the variation of the distribution more 
accurately than the long wavelength. Consequentially, in the com- 
prehensive analysis the experimental 0-X R map became more similar 
to the 0-l R map in Fig. 5(f) than that in Fig. 5(d)) through statistical 
process. The similarity between the experimental and calculated 
graphs can give a validity of 2D simulation at least in this case. 

Discussion 

We acquired the distinguishing 0-l R maps of the PSP-LSP coupling 
in an AgNC@Si0 2 monolayer using 2D simulation to effectively 
process the enormous data. To verify the appropriateness of the 
2D simulation results, representative 3D simulations were con- 
ducted. The angle (<p) between the perpendicular plane (xz-plane) 
to the monolayer and the orthogonally projected incident light on the 
monolayer (xy-plane) was additionally considered in 3D simula- 
tions, as shown in Fig. S6(a) of the Supplementary Information. In 
the rigorous coupled wave analysis (RCWA), the accuracy of numer- 
ical results is proportional to the number of eigenmodes that were 
used for calculation, but it has trade-off relations with computation 
time 19 22 . Therefore, we conducted 3D simulation with lower eigen- 
modes and resolution than those of 2D simulation in the limited 
region of wavelengths from 520 nm to 720 nm and incident angles 
from 55 u to 70°, which is depicted by black lines in Fig. S6(b). As 
results of 3D simulations, the shape and position of two hybrid bands 
on 0-). R map are similar to those of the 2D simulation result as in 



Figs. S6(b)-(d), even though minute differences can be missed due to 
the relatively low number of eigenmodes. Therefore, the similarity 
between 2D and 3D calculations supports the reliability of the 2D 
simulation results. 

The intermediate characteristics of the PSP-LSP coupling in 
QCML lead us to interesting subjects like band engineering and slow 
SPs. In QCML, the electromagnetic fields are changed between the 
intrinsic states of PSPs and LSPs as in Figs. l(g)-(l). This means that 
SPs on the same structure can be minutely tuned by changing the 
excitation conditions. In particular, the lower energy band in 
Fig. 1(d) can be regarded as a new dispersion characteristic modified 
from that of PSPs in a conventional metal surface. In this band, the 
group velocity decreases due to the effects of LSPs. From our numer- 
ical results, the group velocity (v g ) of the SPs in our QCML can be 
lowered down to the less than one fourth times to that of SPs in a 
conventional metal layer. Therefore, it can be potentially applicable 
to design the energy bands in plasmonics or to generate the slow SPs. 

In conclusion, we can hybridize two conspicuous dispersion bands 
with the PSP-LSP coupling in a QCML and concretely demonstrate 
their intermediate behaviors depending on structural parameters. 
Although a precise QCML structure with a constant gap of few 
nanometers cannot be realized yet with neither top-down nor bot- 
tom-up process in current technology, the intermediate characteris- 
tics between PSPs and LSPs on the proposed structure will clearly 
show the distinguishing property if an innovative process for the 
nanofabrication is developed in future. The intermediate character- 
istics have the advantage to simultaneously excite PSPs and LSPs in 
one structure: for examples, because PSP- and LSP-based bio sensors 
having different optimum conditions in sensitivity and resolution for 
target analytes 46,47 can be complementally combined in one chip, the 
complex interactions of individual biomolecules and a bulk cell can 
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Figure 5 | Statistical feedback of spatial distribution of AgNC@Si0 2 s to simulation, (a) Inter-particle distance distribution of AgNC@Si0 2 s in sample. 
The average inter-particle distance is 49.3 nm and standard deviation is 7.47 nm. (b) Randomly generated distribution with normality, (c) Schematic 
diagram of simulation structure and condition for random distribution of inter-particle distances. One period includes one hundred of randomly 
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the same positions on the simulation maps of (d) and (f). 
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be detected by different states of SPs. The control on the energy bands 
and speed of SPs in a QCML can also provide the fundamental 
methods to research nanophotonic signal processing 48 ' 49 . Further- 
more, if it is combined with various flexible or nonlinear materials, 
it will potentially contribute to developing active nanophotonic 
devices because the structural and dielectric properties can be tun- 
able 50 ' 51 to actively utilize the multiple bands and PSP-LSP transition 
in a QCML. 

Methods 

Simulations. For the calculation of electromagnetic fields and reflectance on the 
proposed structure, our own-built rigorous coupled wave analysis (RCWA) tool is 
used 19-22 . The numerical method is firstly to divide overall geometry as an array of 
layers and to calculate the Fourier spectrum of each layer. Then, eigenmodes which 
can exist in those layers are obtained. Finally, coupling coefficients from the incident 
wave to eigenmodes in each layer are calculated by using the extended transfer matrix 
method (ETMM). By summing backwardly the propagating coupling coefficients for 
all the eigenmodes existing in the incident layer, we could obtain reflectance for a 
specific wavelength and incident angle. 9-Xr map is calculated from the repeated 
process of RCWA calculation varying with the condition of incident angle and 
wavelength. 

Preparation of AgNC@Si0 2 s. 50 mg/mL solution of silver nitrate (AgN0 3 , Sigma- 
Aldrich), 25 mg/mL of poly(vinyl pyrrolidone) (PVP, Sigma -Aldrich), and 5 raM 
solution of sodium hydrosulfide (NaHS, Sigma-Aldrich) in ethylene glycol (EG) were 
prepared. 8 mL of EG was heated at 150°C and stirred with a silicon coated magnetic 
bar for 1 h in a 20 mL glass vial. 0.05 mL of the sulfide solution, 1.8 mL of the PVP 
solution and 0.4 mL of the AgN0 3 solution were sequentially injected into the heated 
EG solution. After 3-5 min of reaction time, the solution was quenched with cooling 
water. AgNCs were dispersed in ethanol after three times of rinsing processes of 
centrifugation, dilution, and sonication with acetone and ethanol. To wrap Si0 2 shell 
around AgNC, 25% solution of ammonia and 10% solution of tetraethoxysilane 
(TEOS, Sigma-Aldrich) in ethanol were prepared. 7.6 mL ethanol was stirred in 
20 mL glass vial and then 3 mL of the AgNC solution, 0.2 mL of the ammonia 
solution and 0.6 mL of the TEOS solution were sequentially injected into the ethanol 
in the vial. The mixed solution was stirred at room temperature for 12 h. After 
reaction, excess reactants were washed away in the rinsing process. 

Self-assembled monolayer of AgNC@Si0 2 s. A glass substrate with the hydrophilic 
surface was immersed in a clean Petri dish filled with DI water for AgNC@Si0 2 SAM. 
AgNC@Si0 2 s were dispersed in 1-butanol and then were slowly dropped onto the 
water surface. Once the SAM was formed at the interface of air and water, it was 
settled on the glass substrate by natural drying. 

Optical measurement. Surface morphologies of the samples were measured with a 
field emission scanning electron microscope (FE-SEM; JSM 6700F, JEOL) and a high 
resolution transmission electron microscope (HR-TEM; JEM-3010, JEOL). The 
samples were illuminated from the bottom by white light of halogen lamp with 
changing the incident angle from 40° to 70°. Reflectance spectra were measured using 
a spectrometer (SM240, Spectral Products) at the same angle with that of incident 
light. 
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